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Abstract 

The zero four-momentum and equal mass limits are taken for the bubble diagram 
of scalar fields. It is seen that RTF and ITF are in complete agreement. However 
contributions from this diagram to both retarded and time-ordered functions do depend 
on the order of the limits and can be infinite in some cases. This shows explicitly that the 
relation between the free energy and a derivative expansion of a thermal effective action 
is generally much more complicated that is the case at zero temperature. 



The simplest example in relativistic field theory of a diagram with non-trivial mo- 
mentum dependence is the bubble diagram of figure (|l]). It is therefore illuminating 

Figure 1: The bubble diagram. 

to see exactly what results the standard Feynman rules lead to especially as they turn 
out to be non-trivial in the zero four-momentum, equal mass limit. This limit is also of 
special interest because the zero energy limit is closely tied to the infra-red behaviour 
which plays a vital role in both first and second order phase transitions as, for example, 
the discussion of the cubic term in the free energy for electroweak theory shows (see 
and references therein). 

However the zero energy limit is problematical. In the context of the calculations done 
using the Feynman rules of RTF (Real-Time Formalisms) 0, ^, ||, ||, |^, pathological 
terms of the form [S{K^ - m^)]^^^ appear g |, |, |10|, [r|, fl, One of the 

motivations of this paper is to show that the obvious problems of RTF at zero energy 
do not indicate that RTF is flawed but rather that thermal Green functions at zero 
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energy are intrinsically more complicated than one would expect from experience with 
zero temperature field theory. Further, using the 'simple' scalar bubble diagram as an 



example, it will be clear that both RTF and ITF (Imaginary-Time Formalism) |^, |T5| 
are in complete agreement and show the same difficulties with the zero energy limit. 
It is hoped that this example will clarify the nature of the calculational problems with 
this limit. In doing so we will be extending existing calculations of the bubble diagram, 
rT| , |T2|, |TB|, to cover all possible variations. 



We start by noting the essential properties of two-point Green functions of equilibrium 
field theory as well as their relation to the results of RTF or ITF calculations. This will 
also serve to establish the notation used here. 

In RTF one has a direct handle on time-ordered Green functions. The connected two- 
point function in RTF, 11'^'', is a two by two matrix where by definition the 11 component 
is the real-time time-ordered function, lit 

m^\t, x) = iUtit, x) = Tr{e-^^T0(t, x)0(O, 0)}/Tr{e-^^} (1) 

This is true whatever version of RTF is used (C*-algebra, Thermo Field Dynamics, path 
integral methods, etc.). 

In ITF one is initially calculating the Euclidean time-ordered function, 11/. In terms 
of energy it is only known at imaginary energies which are multiples of 2111/(3. Only by 
making an analytic continuation can one look at the the behaviour of ITF calculations 
in the neighbourhood of zero energy. Using standard boundary conditions at infinite 
energies this can be done [|1^, and the resulting function will be denoted by He. One 



finds that at real energies ITF is calculating the retarded {R) and advanced {A) Green 
functions 

Yic{z = E + ie) = R{E) 

U^{z = E-ie) = A{E) (2) 

where E will be used to indicate real energies. As functions of time, the retarded and 
advanced functions are given by 

tR{t,x) = e(t) Tr{e-^^[0(t,f)</)t(O,O)]}/Tr{e-^^} (3) 
iA{t,x) = -^(-t) Tr{e-^^[0(t,f)0l-(O,O)]}/Tr{e~'3^} (4) 

where only bosonic fields are being considered. 

Standard relations between the time-ordered and retarded functions are simple to 
derive [^] and this then provides a link between the usual results of RTF and ITF calcu- 
lations by using (P and (^. In particular relations can be found for truncated diagrams 



e^^ + 1 

lm{[U-\ET} = -^lm{R{E)}, (5) 

Re{[n-i(E)]ii} = Re{R{E)} (6) 

Im{i?(E = 0)-^} = (7) 



2 



where 



n{E) 



Such relations are derived from the definitions of these full Green functions but they 
hold even when the contribution from only one single Feynman diagram is considered. 
This is not surprising as the relations come directly from the definitions of what is meant 
by equilibrium thermal Green functions (cyclicity of the thermal trace, Kubo-Martin 
Schwinger condition etc.) which must be respected order by order if an approximation is 
to be physically realistic. 

The usual problem of the infra-red divergence in the Bose-Einstein distribution seems 
to hit us straight away and we find for the bosonic case 

\\^^{[Ii-\E)r - R{E)-'] = hm^lm{R{E)-'} (9) 

The divergence in (P) looks very suspicious in view of the RTF divergences we will discuss 
below. However was derived for arbitrary p, mi, 1712 whereas it is the behaviour at 
equal mass and zero four-momentum that will interest us below. Further, as noted above, 
lm{R[E = 0)} = and a more careful analysis shows that 



hra{[U-\E)Y'-RiEr'} = | t/t Tr{e-^^0,(t, x)0t (Q, 0)}/Tr{e-'^^}. (10) 

so that the obvious divergence in has gone and the difference may well be finite. Thus 
does not seem to tell us much about possible divergences. 

The problems which we wish to address here are most evident in the RTF calculations 
where it is well known that singular terms do appear. These are the — m^)]^-^ 

terms which appear when lines in a diagram carry the same four-momentum and are 



associated with particles of the same mass ||, ^, |I0|, |TT], [T^ . This occurs when 

parts of a diagram are self-energy insertions but in this case it is simple to see that when 
the usual RTF sum over internal vertex labelings is performed, such contributions cancel 

i, i ai- 

However such singular terms also arise in diagrams where there are external legs 
carrying zero four-momenta [|10|, |12|. Such diagrams are physically important, for instance 
the free energy or effective potential is also the generating functional of all IPI diagrams 
in which any external legs carry zero four-momenta, though this is not so straightforward 
for the thermal case . 

The simplest example of the [6{K^ — m^)]^-^ singularities comes from the bubble 
self-energy diagram ([l|) where two scalar fields run round the loop. The contribution to 
the time-ordered function, Bt, is easily calculated in RTF and is given by 

- iBt{E,p;mi,m2) = -iB^^ 

= -iBto - iBti - iBt2 

(11) 



where 



ni\h\)5{K' 



m 



(12) 



_ ^2 _|_ 

It is convenient to split it into three terms, Btj being the contribution to the time-ordered 
functions from the terms in the RTF calculation containing j delta functions and Bose- 
Einstein factors. 

Throughout E will represent real external energies with = {E,p)'^ the external 
Minkowskii four-momentum, while K'^ = {k^, k)^ denotes Minkowskii loop momentum. 
The modulus of three- momentum is denoted hj p = \p\,k = \k\. 

The Bt2 contains a product of two delta functions. This is only a problem for the 
RTF calculation if we have set mi = 1712 and E,p = before doing the integrals. Thus 
avoiding any one of these limits is sufficient to give a well defined integral and this is 
what will be done. Doing the energy and angular integrals gives 



Bt2{E,p;mi,m2) 



IGvrp 



E 



k 

dk — n{uji)n{VL2)9{cos 9* 



where 



cos^* 



cos(0) 



1 



2pk 
p.k 
pk 



{P'^ + ml-ml±2Euoi) 



The dispersion relations, cUj, VLi are defined to be 

fc2 + m^, ai = {p + kf + m^. 



1)1^,=^. + (1^2) (13) 

(14) 
(15) 

(16) 



Note that in general the two delta functions can both be non-zero . Thus Bt2 is often 
non-zero and well behaved as ([1^) shows. 

If we keep mi 7^ m2 then as ^ we find 



lim Bt2iE , p; mi ^ m2) = 



;i7) 



as the argument of the two delta functions can not be satisfied simultaneously. However 
setting m = mi = m2 first gives an interesting result. Taking the limit E,p ^ the ratio 
f = liJI can be kept fixed, where 



P_ 
E 
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We then find 



lim lim S+of-E", ^1 7^ ^2! 



—ig mXc 

\p\ 



dx in?{mxcx). (19) 



57r 



T 

¥\ 



1 

'2 ' 



rj.2 



mxc\p\ 



Co(l)-^(l) + log(6)) 



(20) 



where 



Xn 



|f7| 



(21) 



Throughout we shall use ~ to indicate that a high temperature expansion, T ^ E,p, mi, m2 
has been taken but only terms of 0{T) and bigger are given. The integral is generally 
finite but does not have a simple form so for illustrative purposes the high temperature 
limit, T ^ mi, m2 ^ E,p, is given. 

In taking the small E,p, (mi —7712) limits, Bt2, is generally well defined except at the 
E,p = rrii = m2 point itself where Bt2 is infinity or zero. It is also pure imaginary. 
However the residue of the pole depends on how this point is approached in the E,p 
plane. So the delta-squared term of the RTF calculation does refiect genuine infinites in 
the Bt2 term. 

As Bt2 is imaginary, it is useful to look at the remaining temperature dependent 
imaginary parts coming from B^ in case it cancels any of the odd behaviour in Bt2- 



While Bfi does not contain any explicit [6{K'^ 
original motivation for studying Bt2, we find 



m 



')]^ terms, which provided some of the 



Btl{E,p] 1711,1712) 



from which 



9' 



E(2vr)- 



n(fi2) 



2uJi {E±ujiY -nl + ie 



2^2 iE±n2y ~ujf + ie 



(22) 



lm{Bn} 



k 1 



i« E /^^— 0(^(^1) + ^(^2)) ^(cosr 



(23) 



Results are similar to Bt2- First 



lim lm{Bti{E , p; nil ^ 1712)} = 0. 
Then we have in the high temperature limit T 3> mi, m2 ^ E,p 



lim lim \m.{Bti\ 

E,p^O 1711^1712 



-^(^ -1), , 

Stt ^ ' \p\ 



(24) 



(25) 



In fact it is easy to calculate the exact total temperature dependent imaginary part and 
this is 



lim lim Im|i?ti+2l 



^0{v^ - l)^M(3mXc) 
Stt \p\ 



-w 

Stt 
|f7| 



(v^ - 11 



T 



mXc\p\ 2\p\ 



(26) 



(27) 



Both Bti and Bt2 contribute to this divergent result despite the fact that Bn had only 
a single delta function in its integrand. However, away from mi = m2 and E,p ^ 
the result is finite. For instance taking the v = 00 limit in (^Uj) ( i.e. p small but not 
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zero) shows that its divergence only occurs at the zero four-momentum point and is not 
associated with just the zero energy hmit as might suggest. So the delta-squared 
term in Bt2 is therefore flagging a deeper problem, the existence of genuine divergences 
present in the total result for the imaginary part of the time-ordered function, lm{Bt}, 
at the zero four-momentum equal mass point. 

Having established this situation with the imaginary part of the time-ordered result 
as obtained in an RTF calculation, it is interesting to see what is happening with other 
results obtained from the same diagram at the same point. So we will look at the real 
part of the time-ordered function obtained in RTF, at the ITF results from which we 
will extract the advanced and retarded Green functions, and check the raltion between 
them. This will enable us to see if RTF calculations or time-ordered functions are unreli- 
able in this limit and whether ITF calculations or retarded Green functions have similar 
problems. 

For the specific case of the bubble diagram in ITF one initially calculates Bj{27Tin/ (3, p; mi, 
which is only defined at discrete values of the energy, 27rm//?, n ^ Z. The only real value 
of energy which is directly accessible is E = 0. In particular, it should be noted that in 
the ITF calculation setting mi = 1712, E,p = before integrating does gives a unique and 
well defined result. There is no divergence in the integrand as there was with RTF and 
this result will be noted below, It in the other ITF calculations discussed here, the 
mi = m2, E,p = limit is taken after the loop energy integration. 

The ITF function, Bj, can be continued to one defined at general complex energies 
which is denoted by Be- It is found to be 



Bc{z G C,p;mi,m2) = (T = 0) + 



9^ 

2 ± 



n{u!i) 1 ^ n{Q2] 



2uji {z±ujiy-ni 2Q2 {z±n2y-ujf 



{2i 



where z is the complex external energy parameter. The retarded, R{E,p;mi,'m2), and 
advanced, A{E,p]mi,m2), functions correspond to continuations of Be to either side of 



the real energy axis, z — > E±ie |18|. By inspection the ie terms are different from the 
RTF B^^ calculation in which the energy integral has been performed. This allows the 
imaginary part of the ITF calculation to differ from the RTF result. We find for general 
p, nil, ^2 and zero energy that 

Bc{z = + ie,p; 1711,1712) = R{E = 0,p; 1711,1712) = A{E = 0,p; 1711,7112) 

= Bi{z = 0,p;mi,m2) 

= Re{B^\E = 0,p;mi,m2)}. (29) 

The last result is obtained by inspection. These results are consistent with the identities 
(I), and (1). Note that not surprisingly the results obtained using Bj at the zero 
energy point are identical with the case where an analytic continuation is made to general 
complex energies and then — > is taken before any other limit. However, these results 
are questionable if the functions are badly behaved at any point. As we have already 
seen problems when E,p = 0,mi = m2, it remains to be shown that the above relations 
hold at that point. 
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The next stage is to look at the E,p = 0, mi = m2 hmit of the ITF calculation of Be- 
The limit can be taken in various ways and it is easily found that 



lim lim BJE,p;mi,m2) = = 0) + -— ■- — - / dk — n(uj 

l^m2 E,p^O AtT^ dm^ Jo LU 



mi— >m2 E,p^O 



{T = 0)-f-,r du Uiu) (30) 



IGtt m 



(31) 



= lim lim lim BJE,p: mi, mo) (32) 

= lim lim lim BJE,p;mi,m2) 

= lim Bj(z = 0,p;mi,m2) 

= Bj{z = 0,p = 0;mi = m2) (33) 

Note that these B^ results coincide with those from the direct, no analytic continuation 
ITF calculation, i.e. Bj (P^D, however the equal mass and zero three-momentum limits 
are taken. In particular, it agrees with the result obtained when setting mi = m2 and 
E = p = before doing the energy sum which we have denoted as (^) with no limit 
symbols. The only one that causes problems is when mi = m2 is taken first and then we 
take — >■ with or after p — >■ 0. If we are to get a single answer we must take E to zero 
before the energy integral is done or before both of the other limits have been taken. If 



this is not done then more complicated results are obtained and I find [0, ^ 



lim lim Bc(E,p:mi,m2] 

p,E^0,p/E=vmi^m2 



rp' roc Ir 

(T = 0)-f^/ dconiPco) , , , (34) 



9' 1 T 
Von 1 + 7 m 



(35) 



In the case f = oo (i? ^ before p ^ 0) (|33D is recovered. 

This general result can be compared against calculations where either E ox p are kept 
non-zero and not small but in which, apart from this, the usual limit is taken. These 
calculations give 

fp' POO 1^ 

B,{E,p = 0,mi = m2) = (T = 0) - / du; ni(3iu) ^ (36) 



BJE = 0,p,mi = m2) = (T = 0) - -— r— / du n((3uj) lo, 

IbTT^ p Ji 



'p + 2k" 
p-2k 



(37) 



In the zero four-momentum limit the above results do indeed reproduce (p^ ) for v = 
and f = oo respectively. 

We can also look at the imaginary part of Be in the case where mi = m2 is taken first 
(p3). This expression for Be is real for |f | < 1 but for |f | > 1 the integrand has a pole 
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and an imaginary part is generated. The small ie present in the energy, and hence in 
tells us which way round the pole to go and we find 



lim lim ImiBJz = E + irilel.p.niT.mo)} = 9(v^ — l)n(dmXr) — (38) 

Note that this ITF calculation of the bubble diagram's contribution to the retarded 
function in the zero four-momentum and equal mass limit is not divergent unlike the the 
time ordered function. However it is not generally zero as suggests. While it is linear 
in energy near zero energy, as the imaginary part of the retarded function should be P 



there is a competing 1/p factor as p — which keeps the result finite. It is only when 
the f — i> oo limit of (^Sj) is taken do we get zero for the imaginary part of Br and so 
agreement with the spectral representation result (|^). This is not surprising as was 
derived for fixed p and E = which corresponds to the v = oo limit. 

One can then use the spectral function results together with these ITF calculations 
of Be to tell us how the imaginary part of the time-ordered function behaves. It is 
completely consistent with the direct RTF calculation of the time ordered function (^6]). 
In particular the time-ordered bubble diagram is found to be infinite for f > 1. However 
it is not the factor oi 1/E for small energies in that is responsible for the divergence, 
this is canceled by the retarded result being linear in E. Rather it is the 1/p behaviour 
at small p that is leading to the singularities in the time-ordered function. 

By inspection of ( p2D and (|28|) , which differ only in the ie factors, the real parts of the 
time-ordered function calculated using RTF, B^^, and the retarded function calculated in 
ITF from Be, are seen to be equal even in the tricky limit. Thus it is clear that RTF and 
ITF are in complete agreement with each other. If one choose to, one could calculate the 
retarded function in RTF or the time ordered function using ITF for the same diagram, 
and the same results would be obtained. 

On the other hand if we look at a different approaches to the zero four-momentum 
point, V finite, then we see that the retarded function has a non-zero, finite value but the 
old infra red divergence in the Bose-Einstein function ensures that the imaginary part of 
the time-ordered function is blowing up in the region 1 < |f | < oo. It is also zero for 
time-like limits f < 1. 

There is one approach to the E, p, (mi — 1712) limit that has not been considered 
so far and that is taking p ^ first and then taking the remaining limits in either order 
or simultaneously. The result is found to be 

lim lim Bc{E,p; mi, 1712) 
= (T = 0)--^ / dk\—M^) + 



671^ Jo I UJ^ 



(fjSm^ n{ijj)\l + n{uj)\- 



Sll'^ Jm 



duo niuj) { — + 



Cj2 



2 2 ^V + 2^^-3a;W | 
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where 



2(mi — 7712) 
E ' 



(41) 



The g = hmit means that E only after the other two and gives the same answer 
as the V = hmit of (Q). This shows that the relative order of the mi —>■ m2 and 
p —>■ is irrelevant, it only matters how they are taken relative to — > as has be found 
elsewhere. The q = 00 limit is the same as (|30D as it should be. 

Lastly we can take the high temperature limit first which allows us to keep E, p, mi, m2 
arbitrary relative to one another and yet get a closed expression for the integral. The 
leading term is 0{T) and is given by making the replacement n{uj) ^ T/uo (e.g. see 
appendix of ll^). Thus by using |^ 



dy- 



y 



52 



log 



y"^ + 2ay cos(t) + 
y"^ — 2ay cos(t) + 



1 



-TT — nt + 71 arctan 



a,b> 0, <t < Tc 



(a^ - 62) cos(t) 



(a2 + 62) sin(t) + 2ab 



(42) 



we find that 



Bc{E,p,mi,m2] 



B 




+ (1 



2) 



(43) 



P2 4^2 



Comparing this result against the high temperature results already obtained independent 
of (^2]) (only (^0]) was derived directly from (|4^)) we see that there is complete agreement. 

It is therefore clear that both time-ordered and retarded functions are not analytic 
near zero-four momenta when the masses are equal. In particular whenever the zero 
energy limit is last, after or taken with the mi = m2,p limits, various answers can be 
obtained. If the masses are kept different [0, |T6| or three-momentum non-zero till last 
then a unique answer is obtained. The results are summarised in figure @j. Starting at 



Figure 2: The uj integrand of B(. for the different possible orders of the limits. 

the centre the energy integral or sum has been performed. The limits are taken in the 
order indicated as one moves out along a radius. The result obtained for B^ is shown 
at the end of a radius. In figure (^ it is given as the integrand for a remaining uo 
integral where a factor of —g'^I^Ti'^ has also been taken out c.f. (pOD,(|34D,(|39[). The high 
temperature equivalents are given in figure (|^) c.f. (PT|),(P3|),(P(]|). Either RTF or ITF 
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Figure 3: The high temperature results of for the different possible orders of the limits. 

can be used to calculate these results. 

The answer obtained when the E ^ limit is not taken last is of particular note 
because it is also the answer obtained if no analytic continuation is done in ITF, where 
E,p = and mi = m2 are all set at the very start before any integration is done. In this 
case the order of the limits is irrelevant because the integrand is analytic. It is not clear 
why these limits coincide but it means one can find this result by keeping mi ^ m2 or 
p ^ till the end. The results obtained in these cases are the same as the one used in 
free energy calculations, 

The real part of the retarded function is always finite as E,p ^ 0. However the 
imaginary part is not always zero its appearance reflecting the lack of analyticity despite 



the fact that it is linear in the energy, (|38D and |T^. Thus there are peculiarities hiding 
in this limit which are not obvious from the ITF formalism usually used to calculate the 
retarded functions. 

The time-ordered function, B^^ of RTF, is found to have the same real part as the 
retarded and so it is also finite. However the imaginary part shows the effect of the 
problematic S"^ term of RTF which results in the same \p\~^ behaviour as is found in the 
retarded function a.sE,p—>-OifE<p and provided mi = m2. Thus the 5^ reflects a 
genuine divergence in the time-ordered function. It is, however, zero and well behaved 
when mi 7^ 1712 or p 7^ 0. 

The problems with the bubble diagram and its various limits may seem to be mere 
technicalities but they do have real do have real implications for calculations of physical 
quantities. One example of interest is the use of zero four-momentum diagrams to calcu- 
late approximations to the finite temperature effective action A detailed discussion 



given in [^. This might be done in order to obtain some dynamical information such 
as is used in studying the electroweak phase transition (see and references therein). 
Usually a derivative expansion on the effective action is being performed. At zero tem- 
perature the lowest order term is just the effective potential and such an expansion is 
well defined The effective potential is thus generated by the zero-momentum IPI 
diagrams. The various results for the zero four-momentum limit of diagrams at finite 
temperature therefore calls into question the simple derivative expansion of thermal ef- 
fective actions but only when there are equal masses or suitable self-interactions in the 
theory. This means that the coefficients of the expansion depend on the order in which 
the time and space derivative expansions are taken. Thus, the results given here suggest 
that the lowest order term is only necessarily the standard free energy if the time deriva- 
tive expansion was done first i.e. a static limit taken, E ^ first. This does not seem 
to be very relevant to dynamics where one might expect a homogeneous, p — >■ limit to 
be more relevant but the example of the bubble diagram suggests that the lowest order 



term is quite different in this case. The calculations performed here and in |jTT], stress 
the role of equal masses in this problem of zero four-momentum limits. 

Finally note that RTF and ITF give completely consistent results for the bubble 



diagram (c.f. comments in |22, 53, M). Either can be used to calculate any of the 



functions discussed here. While the RTF has obvious singularities in the zero four- 
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momentum equal mass limit, the [5{K'^ — m^)]^-^, the same singularities are hiding in 
ITF calculations of time-ordered functions. Related problems appear in calculations of 
the retarded function in the unexpected non-zero, if finite, imaginary parts. 
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